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. Abstract 

A fission decay of highly excited periodically driven compound nuclei is considered 
in the framework of Langevin approach. We have used residual-time distribution (RTD) 
as the tool for studying of dynamic features in a presence of periodic perturbation. The 
structure of RTD essentially depends on the relation between Kramers decay rate and 
the frequency oj of the periodic perturbation. In particular, intensity of the first peak 
in RTD has a sharp maximum at certain nuclear temperature depending on to. This 
maximum should be considered as first-hand manifestation of stochastic resonance in 



nuclear dynamics. 



1 Introduction 

j> | The atomic nucleus since its discovery has been constantly used for verifying of new physical 
ideas such as tunneling [|IJ, superffuids 0, superconductivity ||, supersymmetry [|J, dynamical 
chaos H . Thus it seems unnatural that one of most recent and intriguing discoveries in nonlinear 
physics-stochasic resonance (SR) (see || for a recent review) up till have not found response of 
the nuclear community. This is particularly odd because there is no doubt that the theory of 
the collective nuclear motion pretending on a consistent description of nuclear dynamics must 
be essentially nonlinear theory. The aim of the present work is to demonstrate the principle 
possibility of observation of SR in nuclear dynamics. As a concrete example we consider a 
process of induced nuclear fission in the presence of weak periodic perturbation. 

SR was introduced nearly 20 years ago to explain the periodicity of the Earth's ice ages 
|7|, U and has found its numerous applications into such diverse fields like physics, chemistry 
and biology (see ||). 

The mechanism of SR can be explained in terms of the motion of a particle in a symmetric 
double-well potential subjected to noise and time periodic forcing. The noise causes incoherent 
transitions between the two wells with a well-known Kramers rate § rt. If we apply a weak 
periodic forcing noise-induced hopping between the potential wells can become synchronized 
with periodic signal. This statistical synchronization takes place at the condition 

r- k l = n/u (1) 



where u is a frequency of periodic forcing. Two prominent feature of SR arises from synchro- 
nization condition ([!]): 

(i) signal-to-noise ratio does not decrease with increasing noise amplitude (as it happens 
in linear system), but attains a maximum at a certain noise strength (optimal noise amplitude 
can be found from ([]]) as is simply connected with it); 

(ii) the residence-time distribution (RTD) demonstrates a series of peaks, centered 
at odd multiples of the half driving period T n = 2(n — |)- with exponentially decreasing 
amplitude. Notice that if a single escape from a local potential well is the event of interest then 
RTD reveals the dynamics of considering system more transparently than the signal-to-noise 
ratio. These signatures of SR are not confined to the special models, but occur in general 
bi- and monostable systems and for different types of noise. 



2 Langevin Equation 

Kramers || was the first to consider nuclear fission as a process of overcoming the potential 
barrier by the Brownian particle. A slow fission degree of freedom (with large collective mass) 
is considered as Brownian particle, and fast nucleon degrees of freedom — as a heat bath. Ade- 
quacy of such description is based on the assumption that the while of equilibrium achievement 
in the system of nucleons degrees of freedom is much less than the characteristic time scale 
of collective motion. The most general way of description of dissipative nuclear dynamics is 



Fokker-Planck equation [K|. However for demonstration of qualitative effects it is convenient 



to use Langevin equation JTTJ that is equivalent to Fokker-Planck equation but is more trans- 
parent. As it has been shown the description based on Langevin equation adequately represents 



nuclear dissipative phenomena such as heavy-ion reactions and fission decay [12], |X3j, [14| and 
possesses a number of advantages over Fokker-Planck description. 

Because we only intend to qualitatively demonstrate SR in nucleus let us consider the 
simplest type of Langevin equation — one- dimensional problem with inertial M and friction 
7 parameters independent on coordinates. Fission coordinate R is considered as a coordinate 
of Brownian particle. The rest degrees of freedom play a role of heat bath being modeled by 
random force £(t). 

The particle motion is described by Langevin equation for canonically conjugate variables 
{P,R} 

dR P 
~dt ~ M 

dP dV , 

_ = -w-^+m (2) 

(3 = 7/M 

£(t) is stochastic force possessing statistical properties of white noise: 

(£(*)) = 0, (£(*)£(*')> = 2L>5(t-t'), D = 1 T (3) 



The nuclear temperature T(MeV)= ^E*/a where E* is an excitation energy and the level 
density parameter a = A/10 (A being a mass number). The (deformation) potential V is given 



as 







_ f 37.46 (R - l) 2 (MeV) for < R < 1.27 

" \ 8.0 - 18.73 (R - 1.8) 2 (MeV) for R > 1.27 



(these are parameters of 205 At nucleus |12||). 

Plausible sources of periodic perturbation are considered below. 
The discretized form of the Langevin equation is given by []13 , |14 



Rn+i — Rn + tPJM 



P n+l = P n {l-(3r) 

Here t n = nr and t](t n ) is a normalized Gaussian-distributed random variable which satisfies 

( V (t)) = 0, (v(tnHtn')) = N5 nn , (6) 

Efficiency of numerical algorithm (|5]) was checked for the following cases: 

(i) V = 0, A = 0, where numerical and analytical results for (P 2 ) and (R 2 ) can be 
compared [jT^|; 

(ii) V ^ 0, A = 0, where numerical and analytical values for Kramers decay rate 
can be compared. According to 0] 



2n 



\Jf3* 2 + 1-/3* 



exp(-AV/T), /T = -^- (7) 



Here cu m i n and a; maa . are the angular frequencies of the potential (J3j) at the potential minimum 
and at the top of barrier respectively, AV is the height of the potential barrier. Numerical values 
of Kramers decay rates for the time bin i is calculated by sampling the number of fission 
events {Nf) i in the i th time bin width At normalized to the number of events N tota [ — J2 (Nf) ■ 

which have not fissioned 

r j = \ Wk (8) 

N total - E {Nf^ At y ' 

Comparison of (0) with asymptotic value of (§) was used for determination of the time 
interval r, which provides saturation for numerical integration (|5]). On the other hand, the 
interval r should be chosen larger than the correlation time of the random process Results 
of numerical calculations are plotted on Fig.l according to (f|) under different number of time 
steps per unit nuclear time h/MeV. One can see that even 20 steps per nuclear time provides 
a sufficient saturation. 



3 Stochastic Resonance in Nuclear Fission 

Now let us proceed to the description of expected effect — manifestation of SR in nuclear fission. 
In the absence of periodic forcing, RTD N(t) has the exponential form (see 0) N(t) exp(— r^t). 
In the presence of the periodic forcing, one observes a series of peaks, centered at odd multiples 
of the half driving period = 2-k/uj. The heights of these peaks decrease exponentially with 
their order number. These peaks are simply explained | T5f| - The best time for the particle to 



escape potential well is when the potential barrier assumes a minimum. A phase of periodic 
perturbation may be chosen in such a special way that the potential barrier V(R) — ARcos(u)t + 
<p) assumes its first minimum at t = 1/2 T^. Thus t = 1/2 T w is a preferred residence time 
interval. Following "good opportunity" to escape occurs in a full period, when potential barrier 
achieves its minimum again. The second peak in the RTD is therefore located at 3/2 T^. 
The location of the other peaks is evident. The peak heights decay exponentially because the 
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Figure 1: The results of numerical calculations for r k under different number 
of time steps per nuclear time h/MeV 



probabilities of the particle to jump over a potential barrier are statistically independent. As 
is shown for symmetric double- well potential |jl6fl , the strength Pi of the first peak at 1/2 
(the area under peak) is a measure of the synchronization between the periodic forcing and 
the switching between the wells. So, if the mean residence time (MRT) of the particle in one 
potential well is much larger than the period of the driving, the particle is not likely to jump over 
the first time the relevant potential barrier assumes its minimum. The RTD exhibits in such a 
case a larger number of peaks where Pi is small. If the MRT is much shorter than the period of 
the driving RTD has already decayed practically to zero before the time 1/2 T w is reached and 
the weight Pi is again small. Optimal synchronization, i.e., maximum Pi, is reached when the 
MRT matches half driving period, i.e., condition (|l]). This resonance condition can be achieved 
by varying the noise intensity D (or u). 

We will show that the same correlation between periodic forcing and escape time takes 
place for a decay of excited states (fission) with a single potential minimum as well. For RTD 
constructing (and following Pi calculation) we use the numerical solutions of Langevin equation 
(H). Let us study evolution of Pi within the temperature interval 1 MeV < T < 6 MeV. 
Corresponding Kramers rates r k and resonance frequency satisfying (fl|) are represented in 
Table 1. Let us fix a frequency of periodic perturbation to = 0.0267MeV/h {T u /2 = 117h/MeV) 
- a resonance frequency at T = 3MeV (see Table 1). On account of the exponential decay 
of peaks heights in RTD (H n ~ exp(— r k T n ), T n = 2(nl/2) tt/uj) , one must observe a series 
of resonance peaks at T < 3MeV. On the other hand, at T > 3MeV (and for the same 
frequency of periodic perturbation) vast majority of nuclei would decay in a while shorter than 
Ti ~ Tu/2. Due to this a sharp maximum of first peak intensity should be observed in the 
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Table 1: 



vicinity of T ~ 3MeV } that is to be interpreted as a manifestation of SR. 

The results of numerical procedure for RTD are presented on Fig.2. Pictures correspond to 
values of Kramers decay rate ( T = 1 — 6 MeV ) under fixed parameters of periodic perturbation 
(A — 1, u — 0.0267). In accordance with expected behavior in the first case (at low ) one can 
distinctly see three peaks located near t = T/2(~ 117.7), 3/2T(~ 353), 5/2T(~ 588), and in the 
second case almost all RTD is concentrated near t = (with width less than T/2). Connected 
with these variations of 1st peak intensity (that represents the measure of the synchronization 
between the periodic forcing and the nuclear temperature and consequently measure of SR) 
are depicted on Fig. 3 for two frequencies of periodic perturbation. Maxima of intensities Pi 
coincide with chosen frequencies of periodic perturbation (see [I]). 
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Figure 2: RTD for T = 1 — 6 MeV 



In conclusion, let us briefly consider the possible sources of periodic perturbation. The 
first possibility is the fissile nucleus as a component of double nuclear system formed, for 
example, in heavy- ion collisions ||17||. In this case, deformational potential will experience 
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Figure 3: Dependence of Pi on T for two different u. 

periodic perturbation similar to tide-waves on the Earth caused by the Moon rotation. In the 
case of asymmetric fission the source of periodic perturbation may be alternating electric field. 
The problem of choice of periodic perturbation would be discussed separately. 
We would like to thank A.Yu. Korchin for valuable discussions. 
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